Controlling the Population Imbalance of a Bose-Einstein Condensate by a 
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Nonlinear Floquet states associated with a symmetry-breaking driving field are exploited to con- 
trol the dynamics of a Bose-Einstein condensate in a double-well potential. The population imbal- 
ance between the two wells is shown to be controllable by slowly tuning system parameters along a 
closed path. The results extend symmetry-breaking-based quantum control to many-body systems 
and extend navigation of linear Floquet states to navigation of nonlinear Floquet states. 
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Symmetry-breaking driving fields have found a variety 
of applications in quantum control, including the gener- 
ation of photo-currents by harmonic mixing [l[, the re- 
alization of basic mechanisms of ratchet transport [3|, 
and the control of chiral processes Q, to name a few. 
To explore how symmetry-breaking driving fields can be 
used to manipulate the dynamics of Bose-Einstein con- 
densates (BEG), it is necessary to extend symmetry- 
breaking-based quantum control scenarios to interacting 
many-body systems. 

Given the robustness inherent in adiabatic quantum 
control, one wonders how a driven BEG system may be 
manipulated by slowly tuning the system parameters, es- 
pecially those of the driving fields. To that end one also 
needs to extend the conventional Floquet states to mean- 
field nonlinear Floquet states jj, la, [H, ll| i ^ new concept 
for understanding periodically driven nonlinear systems. 
Nonlinear Floquet states may display degeneracy as well 
as bifurcation points that are absent in linear systems 
[a, Q- This presents both challenges and opportunities 
for adiabatic quantum control. In particular, the number 
of nonlinear Floquet states can differ from the dimension 
of the Hilbert space, and that the adiabatic following of 
nonlinear systems might necessarily break down due to 
bifurcations [g]. 

In this study we demonstrate that the population im- 
balance of a BEG in a double-well potential [^, llOi llllj 
[12, [l3] can be controlled (including its inversion) by using 
a symmetry-breaking driving field. Specifically, by con- 
sidering a closed path of the system parameters, we show 
that it is possible to move one self-trapped state from one 
well to the other, thus inducing an inversion of the pop- 
ulation imbalance and placing the system in a different 
nonlinear Floquet state. Recognizing that a double-well 
model is relevant to other contexts as well, e.g., a BEG 
occupying two hyperfine levels [Ij] or a BEG in an opti- 
cal lattice that occupies two bands [1^1, this work shall 
motivate further studies of robust control of many-body 
quantum systems with symmetry-breaking driving fields. 

A BEG in a double-well potential can be readily re- 
alized [ill, [l3, [i3- For example, the double- well po- 



tential can be implemented by using two optical lattice 
potentials [16[. The periodic driving force can be cre- 
ated by modulating the relative phase between the two 
optical lattices '13'| (see also [13 )■ Theoretically, under 
the mean-field two-mode treatment [18|, the associated 



driven dynamics can be described as follows: 



(1) 
(2) 



where C denotes the tunneling rate constant, g describes 
the nonlinearity due to the self- interaction of the BEG, 
-^1.2 = |'0i,2p are the relative populations in the left or 
right wells (under the normalization A^i + N2 = 1), 

fit) = A sin(cot) + /2 sm{2Ljt + 9) (3) 

is a zero-bias bichromatic driving force with period T = 
27r/w, and 6 is the relative phase between the lo and 2lu 
fields. For a general value of 9 the bichromatic driving 
force f{t) breaks the time-reversal symmetry and a gener- 
alized parity [2'|. Throughout we set ?i = 1, T = 1 (hence 
all other parameters should be understood as scaled di- 
mensionless variables). Typical values of these system 
parameters were well discussed in, e.g., Refs. [j, [l9| . 
The value of g is assumed to be tunable via the Feshbach 
resonance under an external magnetic field, though the 
magnetic field is not explicitly included in the Hamilto- 
nian. 

For f{t) = 0, Eqs. ([J) and ^ have a left-right per- 
mutation symmetry and reduce to the static double-well 
model for self-trapping [^, [lO, \1M- That is, for a suf- 
ficiently large g, there exist two degenerate eigenstates 
with a nonzero population imbalance S = Ni — A^2- For 
f{t) 7^ 0, in the linear limit the system ([I|)-([l|) pos- 
sesses two Floquet states whose time evolution satisfies 
[MT),MT)] = [</)i(0),(^2(0)]exp(-ie), where 0i,2(i + 
T) — 01,2(0 ^^^ ^ is the quasienergy. The Hnear Floquet 
states can then be continuously extended into the non- 
linear domain (7|, thus defining nonlinear Floquet states. 
Like in the stationary case, Floquet states can also show 
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FIG. 1: (a) Quasienergy of nonlinear Floquet states vs the 
nonlinearity strength g. Points 1 and 2 represent the initial 
and final states in an adiabatic process (arrows indicate the 
direction). The population imbalance of state 3 on another 
branch is almost the opposite of that for state 2. (b) Pop- 
ulations for each well averaged over a driving period as g is 
adiabatically increased. The dashed lines labeled by 1 or 2 
correspond to the states in (a). The thin solid line marks the 
threshold value g^ for the onset of self-trapping for Q = —1.6. 
Insets depict the populations of atoms in the two wells for 
states 1 and 2. Other parameters are uj — 2-k, fi ~ f2 ~ i, 
C= 1. 



population imbalance for large g as a manifestation of 
self-trapping. 

To elucidate the properties of the nonlinear Floquet 
states we may also map the mean-field dynamics of the 
two-mode system onto the dynamics of a periodically 
driven classical system. Introducing the variable if = 
■di — '&i, with ipi = IV'il exp (zi?i) and i/)2 = |V'2| exp (1^2), 
a classical effective Hamiltonian H^s exactly describing 
the dynamics of S and ip is found to be [5| HaS — 
IgS"^ + CVl - 52 cos((p) + 2f{t)S. The nonlinear Flo- 
quet states can then be mapped to the fixed points of 
the Poincare map of HeS 0, Q • 

To examine how the nonlinear Floquet states may be 
adiabatically correlated, we first study their behavior as 
the nonlinear parameter g increases. As we increases g, 
one of the two Floquet states can bifurcate into three non- 
linear Floquet states at a critical value g^- For the special 
cases of = 0, tt, the system conserves a generalized par- 
ity (time reversal together with a left-right permutation), 
we get a pitchfork bifurcation with two degenerate states 
with exactly opposite nonzero population imbalance S. 
When the above symmetry is broken a saddle node [Fig. 
[T] appears instead, where two of the new states have very 
close quasi-energies with almost opposite S] and the third 
state corresponds to a saddle point on the Poincare map 
of i^off • In this latter case the critical nonlinearity value 
for the bifurcation is a function of 0, denoted gc{Q)- The 
arrows in Fig. [l]show that with an adiabatic increase in 
g the system will continuously follow the upper bifurca- 
tion branch corresponding to a nonlinear Floquet state 
with increasing S. 

Taking the parameter g as an example, how slowly the 



system parameters should change in order to achieve nec- 
essary adiabaticity? To answer this question we consider 
cases of linear ramping, where g is made to change lin- 
early in time, between g — gi and g = 52, at a rate of a. 
To get a rough estimate of how small the ramping rate 
a should be to ensure adiabaticity, we use the Landau- 
Zener tunneling formula of Floquet states in the linear 
limit [205. For a representative case such as that shown 
in FiglU we find that a ^ 10"'* suffices for an adiabatic 
process for which the Landau-Zener tunneling probabil- 
ity becomes negligible. To check that we take gi — 1.8 
and g2 = 2.2 and consider one initial state such as the 
state "1" shown in Fig. [1] (preparing such an initial state 
requires the tuning of S and ip [11], followed by a fast 
turn-on of the driving field). After the ramping we let 
the state evolve for some time and then examine the time 
dependence of A^i and iV2- For a < 10^^, the obtained 
state after the ramping process is essentially identical 
with the adiabatic state associated with g = g2 (or state 
"2" in Fig. 1). The periodic oscillations in Ni and A^2 
also exactly match the frequency of the driving field. 

Significantly, if the main interest is in the control of the 
population imbalance S, then the averaged S over many 
cycles of the driving field is more relevant and hence it 
is unnecessary to achieve full adiabaticity. Consider then 
cases with considerably rapid ramping, e.g., a = 10"^ (so 
the ramping only need r^ 10^ field periods). As shown in 
Fig. 2, for such a ramping rate that is two orders of mag- 
nitude larger than before, almost the same time-averaged 
S is obtained. This can be further appreciated using an 
alternative perspective afforded by the classical effective 
Hamiltonian Hcs- As seen from the right panel in Fig. 
2, if a is much larger than 10"^, self-trapping breaks 
down. The classical trajectory after the ramping pro- 
cess is found to be on a chaotic layer around a separatrix 
structure, which corresponds to the state in the lowest 
bifurcation branch. This is understandable because after 
a fast ramping process, the system has not evolved much; 
so if the initial state is the fixed point for g = 1.8 located 
at (^ w —0.3 and 5* w 0, it will necessarily find itself in 
the neighborhood of the saddle point for 5 = 2.2 located 
also at (/? w —0.3 and 5* ~ 0. 

By contrast, if a < 10^^, then the classical trajectory 
remains in the neighborhood of the fixed point. Averag- 
ing out the oscillations around the fixed point we obtain 
an averaged S close to that in the fully adiabatic case. 
Evidently then, manipulation of population imbalance 
may be effectively achieved within a relatively short time 
scale. 

We now ask if the population imbalance S can be in- 
verted. For that purpose the parameter 9 becomes crucial 
because it is responsible for the symmetry breaking here. 
However, for strong self-interaction g > gc, a varying 9 
may have little effect on a large S due to the self-trapping 
effect. Further, a path in the parameter space may not 
allow for adiabatic following if certain loop structure of 
nonlinear Floquet states is encountered (see, for exam- 
ple, case (c) of the lower panel of Fig. 2]). This motivates 
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FIG. 2: Left panel: Evolution of the populations A''i and 
A'^2 after ramping the nonlinearity parameter g with different 
rates a. Right panel: The same dynamics as in the left panel 
plotted on the S-ip plane associated with //off in Eq. (4). 
Parameters are the same as in Fig. [l] The case of a = 
represents the true adiabatic limit. 



us to study how S might be inverted by tuning both 9 
and g. 

Consider then a four-stage control cycle with state 2 
in Fig. 1 being the initial state. In stage I, the value 
of g is ramped from g = 2.2 to g = 1.8, a point where 
degeneracy of nonlinear Floquet states no longer exists 
for any value of 9. In stage II, 9 is ramped from —1.6 
to 1.6. In stage III the value of g is ramped back to 2.2, 
and in stage IV 9 is ramped back to —1.6. Stages III and 
IV are hence the opposite actions of stages I and II, but 
with a different order. As we ramp different parameters, 
we can have different ramping rates, e.g., ag, ag. For 
simplicity, we set the two ramping rates to be the same. 

Figure [3] compares the time dependence of S associ- 
ated with state 2 (bottom curve in the left panel) and 
that associated with the final state (upper curve with os- 
cillations similar to state 2) after the four-stage ramping 
cycle, with the linear ramp rate in both g and 9 given 
by a = 10"''. Clearly, the resultant state after following 
the parameter cycle inverts S. Moreover, its dynamics 
coincides with that of state 3 (see Fig. 1). This indi- 
cates that after the four-stage control cycle, a nonlinear 
Floquet state different from the initial one is reached. 
Next we show in Fig. [3] the quasienergy along the con- 
trol path. At first glance, the quasienergy path seems 
to form a closed loop. However, on a closer inspection 
of the last stage of the cycle (see panel (b) on the right 
of Fig. 3), it is found that the quasienergy of the final 
state (square) coincides with that of state 3 (see Fig. 1) 
and hence differs from that of the initial state (circle). 
This again confirms that the system is placed on a dif- 
ferent nonlinear Floquet state after following the above 
ramping cycle. 

To shed more light on the results we calculate the 
quasienergy surfaces as a function of 9 and g [Fig. 0]. 
For the sake of clarity, we only show the surface around 
the region where one quasienergy surface bifurcates into 
three surfaces. Above a critical strength of g, two quasi- 
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FIG. 3: Left panel: Population imbalance S of the initial 
state 2 (bottom curve) and the final state (upper curves) after 
a four-stage ramping cycle, a = 10"'* for the upper solid 
curve and a = 10~^ for the upper dotted-dashed curve. Inset: 
the four-stage cycle with the dot being the starting point. 
Right panel: (a) Quasienergy along the four-stage adiabatic 
cycle, (b) Enlargement of the panel (a). The initial and 
final points are depicted by a circle and square, respectively. 
Arrows indicate the direction of the ramping process. 



energy surfaces can cross each other, forming a degener- 
acy line. Because this degeneracy line is along the 9 = 
direction (consistent with the fact that f{t) with 9 = 
conserves a generalized parity discussed above) , the crit- 
ical value g for the surface crossing can be identified as 
gc{9 = 0). Note also that the degeneracy line does not 
present a problem for the system to continuously follow 
the ramping process, because (i) S can only change con- 
tinuously and (ii) the involved degenerate states however 
have opposite values of S. 

Using the quasi-energy surfaces and its projection at 
different g values, we can now visualize the dynamics 
associated with the above four-stage cycle. In the first 
stage, the state slides down the quasi-energy surface 
shown in Fig. 4, and at the same time three nonlin- 
ear Floquet states are merged into one. At the end of 
the first stage, the self-trapping phenomenon has disap- 
peared and this makes it possible to reverse S by tun- 
ing 9. The second step does exactly this by allowing the 
symmetry-breaking driving force to change S from a very 
small, negative value to a very small, but positive value. 
In the third stage, the g value returns to its original value, 
the system climbs up the quasi-energy surface, and three 
nonlinear Floquet states emerge again as a result of bifur- 
cation. During this stage the inverted S is then magnified 
in the same mechanism as shown in Fig. 2. In the last 
stage when 9 is ramped back, the system can maintain its 
population imbalance due to self-trapping. In the mean- 
time the system crosses the degeneracy line, returns all 
system parameters to their initial values, and finally ends 
up in a different Floquet state. This analysis also makes 
it clear that by considering a reversed closed path one 




FIG. 4: Upper panel: Quasienergy surface vs the relative 
phase parameter 6 and the nonlinear parameter g. Lower 
panel: Quasienergy e vs 6 for various fixed values of (?.(a) 
g = 2; (b) g = 1.98; (c) g = 1.97; (d) g = 1.95. The ini- 
tial state is depicted by a black dot in panel (a). In the first 
stage of the control cycle, g decreases and the quasienergy 
state (black dot) changes as shown in the a-b-c-d sequence 
of the four panels. Arrows in panel (d) indicates the evo- 
lution of the quasi-energy in the second stage. This yields 
a state with almost opposite population imbalance depicted 
by the square. In the third stage, in the order of d-c-b-a, 
the quasienergy (square) increases. The evolution in the final 
stage is indicated by the arrows in panel (a). 



can also change S from positive to negative. 

One may regard the above "anholonomy-like" behavior 
as a consequence of adiabatic following during the entire 
ramping process, insofar as the system always follows the 
control cycle and never populates two nonlinear Floquet 
states. However, due to the nonlinearity-induced degen- 
eracy line along = 0, it is better to regard the evolution 
during stage IV as a diabatic process. That is, during 
stage IV the change of the parameters is always "too 
fast" as compared with the zero spacing at the degener- 
acy point. Indeed, in a fully quantum treatment with- 



out the mean-field approximation, the quantum levels, 
though nondegenerate, will be highly clustered around 
the mean- field degeneracy point. Therefore in a fully 
quantum picture adiabatic following breaks down during 
stage IV. In either perspective, the net result is the same: 
the system finds itself in a different state after following 
a ramping cycle in the parameter space. 

We next examine whether we can achieve an inver- 
sion of S if the system parameters are tuned at a much 
faster rate? Significantly, as already suggested by the 
results in Fig. 2, the adiabatic following approach here 
still works even if we ranrp the parameters g and 9 at 
a considerable rate. Indeed, the result for a ramping 
rate two orders of magnitude larger (a = 10~^) is also 
shown in Fig. 3 (dotted-dashed curve on the left panel). 
In this case, as indicated by the beating patterns of the 
dynamics, the final state is not on one single Floquet 
state, but upon a time average we still achieve an in- 
verted S close to that obtained from the previous slow 
ramping case. We have also checked many other clock- 
wise (counter-clockwise) closed paths in the g — 9 plane 
of Fig. 4, finding that they can yield similar control from 
state 2 to state 3 (from state 3 to state 2) if they enclose 
the critical point g — gc{9 = 0); 6* = 0. Thus, we conclude 
that, given the topology of quasi-energy surfaces, espe- 
cially the degeneracy and bifurcation points, navigation 
of nonlinear Floquet states provides a useful control sce- 
nario. This also extends an early control scenario based 
on the navigation of linear Floquet states j2l|] . Applying 
our approach to the cases with nonlinear Floquet states 
localized in momentum [7|, it should be possible to in- 
vert ratchet transport in many-body systems without a 
biased field. 

In conclusion, by demonstrating the control of the pop- 
ulation imbalance associated with a BEG in a double- well 
potential, we have shown that the navigation of nonlinear 
Floquet states associated with symmetry-breaking driv- 
ing fields can offer effective and robust control over many- 
body systems on the mean-field level. 
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